AD-751  940 


IDENTIFICATION  OF  EARTHQUAKES  AND  UNDER 
GROUND  EXPLOSIONS 

i 

Eugene  Herrin 

Southern  Methodist  University 


s - 

Prepared  for: 

Air  Force  Office  of  Scientific  Research 
Advanced  Research  Projects  Agency 

•  i 

October  1972 


I 


I 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


AD  751940 


24  OCT  1372 


TECHNICAL  REPORT 
to  the 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 

from 

Eugene  Herrin 

Dallas  Geophysical  Laboratory 
Southern  Methodist  University 
Dallas,  Texas  75222 


University 


ARPA  Order  Number:  1827 

Program  Code:  1F10 

Name  of  Contractor:  Southern  Methodist 

Effective  Date  of  Grant:  July  1,  1971 

Grant  Expiration  Date:  June  30,  1972 

Amount  cf  Grant  Dollars:  $90,000 

Grant  Number:  71-2133 

Principal  Investigator  and  Phone  Number: 

Program  Manager  and  Phone  Number : 


Title  of  Work:  Identification  of 

University  Account  Number:  80-46 


Eugene  Herrin 
#  214  692-2760 
Truman  F.  Cook,  Director  of  Research 
Administration 
#  214  692-2031 

Earthquakes  and  Underground  Explosions 


Sponsored  by 

Advanced  Research  Projects  Aeencv 
ARPA  Order  No.  1827  X 


Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 


U  S  Department  of  Commerce 
Springfield  VA  22151 


Approved  for  public  release; 
distribution  unlimited. 


1 


UNCLASSIFIED  _ 


■vnmly  Cl  ■SHincAtipn  ■  ' 

- -  DOCUMENT  CONTROL  DATA  •  R  8.  D 

-S,cu,.„  Casa, „ea, ion  o,  hodv  o'  ahslrar,  and  indesin,  -n  «*,"*  *«„  ^entered 

)'.  orTgina  ting  activity  /Carroratr  author)  ,TMPT  aqqtPTFD 

Southern  Methodist  University  _ UNCLAbSltlbu - 

Dallas  Geophysical  Laboratory  **■  GF,ouP 

Dallas  TX  75222  _ _ _ 

I'  ^Identification  of  Earthquakes  and  Underground  Explosions 


4.  DESCRIPTIVE  NOTES  fDpc  o  t  report  and  inclusive  dales) 

Scientific . Interim 

5.  au  THORtSi  (t  irsf  name,  middle  Initial,  last  name) 

Eugene  Herrin 


le  RF.POF(i  DATE 


i  October  1972 

««.  CCNlMCf  OR  GRANT  NO. 

AFOSR-7 1-2133 

b  PROJEC T  NO 

AO  1827 

C. 

62701D 


7a  TOT  AL.  NO.  OF  PAGES  [7b-  NO  OF  REPS 

>*r  3 )  [_  o 

Ufl.  ORIGINATOR'S  REPORT  NUMBERS' 

80-46 


~^ToTh  F  “rep O R  T  n 0 151  (Any  other  numhera~lhal  may  bo  assigned  ' 

,hU,'PlwSR-TR-r2-2l6  7  ! 


10.  DISTRIBUTION  STATEMENT 


11  SUPPLL.MNTAmV  notes 


TECH,  OTHER 


rTH  ABSTF»»C  I 


Approved  for  public  release;  I 

distribution  unlimited.  _ _} 

- ““ - -  12.  SPONSORING  MILITARY  ACTIVITY  1 

Air  Force  Office  of  Scientific  Research/NPG 
1400  Wilson  Boulevard 

Arlington,  VA  22209  _ _ 


During  the  period  January  1970  through  August  1970  long  period  seismographs 
vS  operated  at  the  Grand  Saline  Salt  Mine  in  North-Central  Texas  This  report 
concentrates  on  the  f clicking  tasks  as  outlined  in  the  original  statement  0 

a  Obtain  digital  data,  in  a  standard  format,  of  long  period  signals  and 
noise  and  to sMre  and  make  available  copies  o£  this  data  to  other  users  on  requestl 

b.  Investigate  methods  of  filtering  which  can  be  used  co  improve  signal  to 
noise  ratios. 

c  Study  regional  variations  in  frequency  dependent  absorption  of  long 
period  surface  waves.  Investigate  the  potential  effect  of  such  regional  variations, 
on  the  estimation  of  surface,  wave  magnitude  and  on  the  estimation  of  the  spect  a 
distribution  of  surface  wave  energy. 
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INTRODUCTION 


During  the  period  January  1970  through  August  1970,  long  period  seis¬ 
mographs  were  operated  at  the  Grand  Saline  Salt  Mine  in  North-Central  Texas 
under  AFOSR  contract  number  71-2133.  Sponsored  by  ARPA  Order  Wo.  1827.  This 
report  concentrates  on  the  following  tasks  as  outlined  in  the  original  state¬ 
ment  of  work: 

1.  Obtain  digital  data,  in  a  standard  format,  of  long 
period  signals  and  noise  and  to  store  and  make  avail¬ 
able  copies  of  this  data  to  other  users  on  request. 

2.  Investigate  methods  of  filtering  which  can  be  used  to 
improve  signal  to  noise  ratios. 

3.  Study  regional  variations  in  frequency  dependent  ab¬ 
sorption  of  long  period  surface  waves.  Investigate 
the  potential  effect  of  such  regional  variations  on 
the  estimation  of  surface  wave  magnitude  and  on  the 
estimation  of  the  spectral  distribution  of  surface 
wave  energy. 

The  seismographs  were  of  the  advanced  long  period  type  with  a  peak  response 
at  50  seconds.  One  three  component  system  was  placed  at  the  surface  and  one 
at  a  depth  of  approximately  180  m  in  the  mine.  The  seismometers  were  sealed  in 
such  a  manner  that  any  "noise"  signal  could  definitely  be  shown  to  be  due  to 
ground  motion  and  not  to  the  buoyant  response  of  the  seismometer. 


3 


r 


DATA  ANALYSIS 


Time  series  of  500  earthquakes  selected  from  NOAA  PDE  Listings,  were 
stripped  from  the  Grand  Saline  field  tapes  and  plotted.  These  events  are 
catalogued  and  stored  on  a  master  library  tape,  readily  available  for  such 
digital  processing  as  filtering  or  spectral  analysis. 

Analysis  of  the  Grand  Saline  data  has  shown  that  significant  amounts  of 
long  period  Rayleigh  wave  energy  were  recorded  from  earthquakes  all  over  the 
world,  with  the  exception  of  those  propagated  along  paths  through  the  Basin 
and  Range  Province  of  Ncrth  America.  Examples  are  shown  in  Figures  1  through 
15.  Figures  1,  2,  and  3  show  in  each  case  signals  from  earthquakes  that  were 
recorded  at  our  Grand  Saline  site  and  the  same  signals  recorded  at  a  site  at 
Queen  Creek,  Arizona.  Gains  are  equalized  at  30  seconds  in  Figure  1  and  at 
50  seconds  in  Figures  2  and  3.  Figure  1  displays  a  small  earthquake  from  the 
South  Pacific.  Note  that  with  the  exception  of  the  superimposed  long  period 
noise  these  signals  are  very  similar.  Both  stations  lie  approximately  on  the 
same  great  circle  arc  through  the  epicenter;  thus,  the  travel  path  to  both 
stations  is  nearly  identical.  Figure  2  shows  an  earthquake  from  the  mid- 
Atlantic,  again  approximately  on  a  great  circle  arc  through  the  recording  sta¬ 
tions.  In  this  case,  however,  the  waves  travelled  from  the  East.  It  is  clearly 
evident  that  the  record  made  at  Grand  Saline  shows  greater  development  of  the 
long  period  surface  waves.  Figure  3  shows  an  earthquake  north  of  Greenland 
recorded  at  the  Texas  and  Arizona  sites.  The  waves  recorded  at  Grand  Saline 
came  essentially  through  the  mid-continent  region  of  the  United  States  while 
the  waves  that  came  to  the  station  in  Arizona  crossed  the  mountainous  region  of 
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Figure 


the  western  United  States.  Note  that  the  longer  period  waves  are  much  smaller 
in  the  Arizona  recording.  Figures  4  through  15  are  earthquakes  which  all  have 
a  travel  path  from  epicenter  through  the  Mid-Continent  region  of  the  United 
States  to  Grand  Saline.  These  figures  show  the  very  long  period  waves  recorded 
at  this  station. 

Fig.  4  mb  4.4,  Ms  4.1  normal  depth,  distance  58°,  source  north  of 

Svalbard,  on  an  azimuth  from  Novaya  Zemlya  to  North-central  Texas. 

The  initial  Rayleigh  wave  arrival  has  a  period  of  about  65  seconds. 
Note  the  dispersed  wave  train  from  the  group  velocity  maximum  at 
about  70  seconds  to  the  minimum  (Airy  phase)  at  about  15  seconds. ^ 

Fig.  5  mb  4.3,  Ms  3.5  normal  depth,  distance  55°,  Iceland,  an  optimum 

path  (approximately  the  same  azimuth  as  that  from  Semipalatinsk  to 
North-central  Texas) .  First  arriving  Rayleigh  wave  has  a  period  of 
greater  than  60  seconds. 

Fig.  6  mb  5.2,  Ms  5.4,  depth  13  km.  distance  121°,  Sinkiang,  note  the  pulse¬ 
like  beginning  at  about  70  sec.  This  event  has  an  unusually  large 
Rayleigh  wave. 

Fig.  7  Energy  density  spectrum ,  of  the  Sinkiang  event  shown  in  Fig.  6. 

This  and  the  following  spectrum  are  not  corrected  for  instrument 
response.  Note  the  significant  energy  at  about  60-70  sec.  and  the 
spectral  notch  at  about  20  sec. 


1.  All  hypocentral  information  and  body  wave  magnitudes  in  this  and  following 
discussion  are  from  NOAA.  Ms  was  calculated  from  Grand  Saline  records. 
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Fig.  8  mb  4.5,  Ms  3.4,  depth  43km,  distance  97°,  USSR -Mongolia  border, 
data  was  high-pass  filtered  with  corner  at  50  sec.  to  enhance 
detection.  Starting  period  is  greater  than  40  sec. 

pig,  9  mb  4.4,  Ms  3.4,  depth  normal,  distance  97°,  Caucasus  high-pass 
filtered  as  for  Fig.  8.  This  event  is  very  near  threshold  for 
single  LP  vertical  instrument  at  Grand  Saline.  First  periods  ob¬ 
served  are  30  to  40  sec. 

Fig. 10  Energy  density  spectrum,  mb  5.0,  Ms  4.2,  Caucasus.  This  figure 

shows  the  spectrum  for  a  larger  event  from  the  same  source  area  as 
for  Fig.  9.  The  spectral  content,  with  significant  energy  to  per¬ 
iods  of  60  sec.  seems  to  be  typical  of  this  source  area. 

Fig. 11  mb  4.5,  Ms  3.5,  depth  25  km,  distance  85°,  Albania.  The  record  is 
underlined  where  the  Rayleigh  wave  was  predicted  to  arrive.  The 
first  observed  period  is  about  50  sec. 

Fig. 12  mb  4.9,  Ms  3.9,  depth  normal,  distance  121°,  Yunnan,  China ,  start¬ 
ing  period  is  greater  than  40  sec. 

Fig. 13  mb  4.9,  Ms  4.5,  depth  normal,  distance  86°,  Lake  Baikal  region.  This 
event  is  mixed  with  the  tail  of  a  Rayleigh  wave  which  arrived  50  min. 
earlier  from  Mongolia  (mb  5.0). 

Fig, 14  mb  4.7,  Ms  3.7,  depth  normal,  distance  70°,  Kamchatka.  The  surface 
wave  travelled  a  rather  poor  path  through  the  Aleutian  arc  and 
northwestern  Canada.  The  starting  period  is  about  60  sec. 

Fig. 15  mb  4.9,  Ms  3.5,  focal  depth  near  100  km,  distance  114°,  Taiwan  re¬ 
gion.  The  first  arriving  Rayleigh  waves  have  periods  of  70  to  75 
seconds.  Note  the  pulse-like  appearance  of  the  Rayleigh  wave  train 
and  the  relatively  low  amplitude  of  the  shorter  period  components 
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(as  would  be  expected  from  an  earthquake  with  this  depth  of  focus) . 
Detection  of  this  event  is  definitely  aided  by  the  enhanced  long 
period  response  of  the  Grnnd  Sflline  seismic  system* 

In  each  case  these  figures  demonstrate  that  long  period  Rayleigh  waves  are  very 
efficiently  propagated  to  northcentral  Texas  along  pr ths  from  northeast,  north 
and  northwest  of  the  station.  We  feel  these  studies,  using  a  single  seismograph, 
have  established  that  north  Texas  lies  at  the  southern  end  of  a  waveguide  extend¬ 
ing  from  the  southern  borders  of  the  USSR  and  northern  China  across  the  Arctic 
Ocean  and  the  Canadian  Shield.  In  the  United  States  the  waveguide  is  bounded  on 
the  west  approximately  by  a  line  through  Carlsbad  and  the  Black  Hills  and  on  the 
east  by  the  Atlantic  and  Gulf  Coastal  Plains.  This  waveguide  is  characterized 
by  relatively  low  temperatures  in  the  lower  crust  and  upper  mantle  which  results 
in  an  exceptionally  good  propagation  path  for  Rayleigh  waves  with  periods  from 
30  to  60  sec.  An  approximate  detection  threshold  of  Ms  3.5  has  been  established 
for  a  single  station  in  North  Central  Texas  for  events  in  the  waveguide  at  dis¬ 
tances  as  great  as  100°.  Because  most  of  the  detections  were  made  at  periods 
greater  than  AO  sec.,  this  magnitude  is  equivalent  to  about  mb  A. 5  for  earth¬ 
quakes  less  than  50  km  in  depth. 

We  have  observed  the  long  period  "noise  hole"  or  "window"  which  has  been 
described  by  using  data  recorded  at  Ogdensburg,  New  Jersey,  Queen  Creek,  Arizona, 
and  Las  Cruces,  New  Mexico.  This  spectra'  hole  is  shown  in  Figure  16  for  both 
the  seismometer  in  the  mine  and  the  surface  seismometer  (during  a  near  calm  per¬ 
iod)  .  It  is  clear  that  any  long  period  seismic  array  should  take  advantage  of 
this  noise  minimum  for  detection  purposes. 
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COHERENCE  SQUARED 


To  establish  an  earthquake  detection  capability  for  the  Grand  Saline  site, 
records  from  the  vertical  long  period  seismograph  were  examined  for  a  30-day 
period  during  May  and  June,  1970.  All  discernible  signals  were  listed  by 
arrival  time  and  period  of  initial  Rayleigh  wave  motion.  It  uas  not  possible 
to  determine  the  location  and  magnitude  of  most  of  these  events  because  they 
mere  not  reported  by  NOAA.  Table  1  shows  the  breakdown  of  the  period’ measure¬ 
ments.  Me  sould  expect  that  close  to  50%  of  those  events  having  a  magnitude 
mb  4.0  and  above  would  be  recorded  at  Grand  Saline  and  of  these,  80%  would  be 
detected  having  initial  Rayleigh  motion  at  greater  than  40  sec.  period. 

Table  2  displays  the  results  of  an  examination  of  Grand  Saline  data  speci¬ 
fically  conducted  to  detect  events  from  the  Slno-Sovlet  Region.  These  data 
indicate  that  the  detection  threshold,  with  relatively  high  confidence,  for  the 
Grand  Saline  station  for  Rayleigh  waves  from  Slno-Sovlet  events  is  approximately 

Ms  3.5  (mb  4.5  for  normal  earthquakes). 

Me  conclude  from  Table  2  and  from  the  previously  shown  seismograms  that 

there  is  significant  Rayleigh  wave  energy  up  to  periods  of  60  sec.  propagated 
along  paths  from  the  USSR  and  Chin,  to  north-ceutral  Texas  and  that  it  is  possible 
to  exploit  the  "noise  window"  centered  at  about  35  seconds  to  improve  Rayleigh- 
wave  detection.  Since  most  events  along  this  waveguide  have  measurable  energy 
at  periods  greater  than  40  seconds  the  Rayleigh  wave  detection  threshold  would 
remain  about  the  same  for  earthquakes  as  deep  as  50  to  60  km.  Surface  wave 
magnitudes  (M,)  reported  here  were  measured  at  or  near  periods  of  20  seconds. 
Magnitude  measurements  from  Grand  Saline  records  based  on  longer  period  Rayleigh 
waves  may  be  more  stable  than  the  conventional  Ms  measurements  because  these 
longer  periods  are  less  affected  by  depth  of  focus  and  small-scale  path  variations. 
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Table  1 


Detection  Study 
30  Days 

During  May  and  June,  1970 
Single  Vertical  LP  -  Grand  Saline,  Texas,  Mine 


Longest  Period 
Detected 


60-70  sec. 
50-60 
40-50 
30-40 
20-30  sec. 


Percentage  of 
Events  Detected 


157. 

387. 

317. 

87. 

87. 


Total  number  of  events  detected  -  380 
Avg.  number  of  events /day  -  12.7 
Predicted  number  of  events/year  -  4623 
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1st  Period  NOAA  G.S.  Signal-to- 

NOAA  Distance  Detected  noise 
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EFFECTS  OF  DEPTH  OF  FOCUS 

I  |  ! 

Theory  based  on  seismic  sources  in  a  layered  half-space  gives  predictions 

,  i 

of  the  effect  of  source  depth  on  the  spectra  of  resulting  Rayleigh  wave?.  In 

'  i 

general,  we  expect  the  spectral  energy  to  shift,  in  a  relative  sense ,• toward 

,  i 

the  longer  periods  as  the  source  depth  is  increased.  Details  of  the  spectral 
shape  depend  in  a  complicated  way  upon  the  nature  of  the  source,  propagation 
azimuth  relative  to  source  symmetry,  the  physical  properties  of  the  layering, 
and  the  exact  depth  of  focus. 

Our  studies  of  Rayleigh  wave  spectra  obtained,  from  the  Grand  Saline  data 

!  I 

show  that  for  many  earthquakes  in  the  Pacific  island  arc  regions  (e.g,,  Tonga, 
Islands) ,  there  is  no  discernable  shift  in  spectral  energy  with  increasing  depth 
of  focus.  Shorter  per iod  Rayleigh  wave  energy  is  generated  jhst  as  efficiently 
by  sources  at  250  km  as  by  those  20  km  deep.  The  theory*  based  on  a  layered 
half-space,  is  clearly  not  applicable  to  these  island  arc  regions,  which  are 

the  source  areas  of  about  90%  of  all  earthquakes.  A  more  realistic  Earth  model,1 

, 

in  which  a  downgoing  slab  of  lithosphere  is  considered,  greatly  compl|icates  the 

theoretical  problem.  Numerical  solutions  to  the  laterally-heterogeneous  pro- 

,  i 

blem  are  possible  in  principle,  but  would  be  difficult  and  time  consuming  in 
practice. 

As  a  first  approach  to  this  problem  we  intend  to  design  and  construct  two- 

i 

dimensional  analog  models  of  an  assumed  island  arc  source  region  with  a  down- 

i 

going  slab.  The  models  will  be  scaled  such  that  Rayleigh  wave  energy  will  be  in 

I  ! 

the  range  20  to  200  kilocycles.  Sampling  techniques  will  be  employed  to  obtain 
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srismic  records  and  to  convert  the  signals  into  digital  form  for  storage  on 
magnetic  tape.  The  model  signals  can  then  be  processed  In  exactly  the  same 
way  as  the  data  from  the  long  period  seismographs. 
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Canton  Experiment 


The  succesr  of  the  Grand  Saline  experiments  led  to  the  installation  of 
a  long  period  triaxial  seismometer  in  a  300  m  borehole  a"  a  site  near  Canton 
This  seismograph  system  was  again  operated  in  conjunction  with  a  3-component 
system  at  the  surface  and  with  tripartrite  5  km  microbarograph  array.  The 
borehole  enabled  the  seismometer  to  be  clamped  at  different  depths.  This  site 
was  in  operation  for  the  period  February  1971  through  May  1972. 

During  this  operating  period  the  TC-200  Data  Acquisition  System  was  record¬ 
ing  for  7870  of  the  time.  Most  down  time  can  be  accounted  to  power  failures  with 
the  rest  under  the  categories  of  maintenance  and  tape  changes. 

Preliminary  analysis  of  approximately  100  selected  events  recorded  at  the 
Canton  site  bear  out  the  findings  in  the  Data  Analysis  Section  of  this  report. 

We  are  currently  in  the  process  of  selecting  events  and  updating  our  event 
library  tapes  as  with  the  Grand  Saline  data.  Fourier  analyzing  will  be  accomp¬ 
lished  for  all  Rayleigh  wave  signals  which  have  a  sufficiently  large  signal-to- 
noise  ratio..  The  spectral  information  is  also  stored  on  a  library  tape  for 
subsequent  processing. 

Digital  recording  of  seismic  data  and  subsequent  digital  processing  tech¬ 
niques  developed  at  Southern  Methodist  University  enabled  us  to  support  the 
LRSM  group  at  Canton  (Teledyne  Geotech)  in  their  studies  of  noise  sources  of  the 
down-hole,  triaxial  long  period  seismometer.  These  same  techniques  along  with 
microbarograph  data  from  our  existing  infrasonic  research  program  supported 
Dr.  G.  (Mike)  Sorrells  (Teledyne  Geotech)  in  his  continuing  measurements  of  long 
period  atmospheric  and  seismic  background  noise,  in  order  to  investigate  the 
causes  of  the  seismic  background  noise  and  perhaps  to  develop  methods  to  min¬ 
imize  that  noise. 


Future  Studies 


Data  from  the  Grand  Saline  and  Canton  experiments  support  the  existence 
of  a  polar  wave  guide  for  long  period  Rayleigh  waves.  A  third  experiment, 
designed  to  collect  still  more  long  period  data,  is  being  planned  for  t  site 
about  35  miles  north  of  Dallas.  At  this  site  we  will  investigate  the  poss¬ 
ible  use  of  information  from  a  multi-element  microbarograph  array  to  pre¬ 
dict  and  possibly  reduce  long  period  seismic  background  noise. 
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